Collective Spin Fluctuation Mode and Raman Scattering in Superconducting Cuprates 
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Although the low frequency electronic Raman response in the superconducting state of the cuprates 
can be largely understood in terms of a d-wave energy gap, a long standing problem has been an 
explanation for the spectra observed in the A\ g polarization orientations. We present calculations 
which suggest that the peak position of the observed Ai g spectra is due to a collective spin fluctu- 
ation mode. 

PACS numbers: 74.72.-h, 78.30.-j 71.10.-w 



X 



Electronic Raman scattering has proven to be a use- 
ful tool in exploring the superconducting state in the 
cuprate materials. The possibility of probing selectively 
electronic excitations in different regions of the Brillouin 
zone (BZ) by the choice of polarization geometries has 
allowed to explore the superconducting gap anisotropy. 
The successful explanation of the Raman data in B\ g and 
B 2g scattering geometries jl],|| has provided one piece of 
evidence for the by now widely accepted d x 2_ y v pairing 
symmetry in hole-doped cuprate superconductors Q . In 
the context of impurity effects as a testing ground for 
unconventional superconductivity, the observed uj 3 to u> 
crossover in the low frequency B\ g Raman response fits 
consistently with the power law crossovers at low temper- 
atures in the NMR relaxation rate and in the magnetic 
penetration depth. An even quantitatively consistent pic- 
ture of electronic Raman scattering and infrared conduc- 
tivity was achieved when the T-matrix approximation in 
the "dirty" d-wave scenario is extended to include a spa- 
tial extension of the impurity potential Q . 

However up to now the discrepancy between Raman 
data in Au, and B\ g ,Bi g geometries has remained un- 
resolved Hlg]. Previous results for the Ai s scattering 
geometry were found to be very sensitive to changes in 
the Raman vertex function 7(k) |^] making a compre- 
hensive explanation difficult for the experimental data in 
different cuprate materials. 

In this paper we calculate the Raman response of a 
d x 2_ y 2 superconductor including the contribution from 
a collective spin fluctuation (SF) mode which is identi- 
fied with the (tt, tt, 7r)-resonance (in short 7r-resonance) 
near lor ~ 41 meV observed by inelastic neutron scatter- 
ing (INS) on bilayer cuprates Our results suggest 
that the A\ g peak position is largely controlled by the 
strength and frequency of the 7r-resonance mode which 
on the other hand does not affect the Raman response 
in the B\ g and B 2g channels. The inclusion of the col- 
lective SF mode allows for a simultaneous fit of the Ra- 
man data in all channels in optimally doped materials. 
Furthermore we find that the inclusion of the SF term 
significantly reduces the sensitivity to the special choice 
of the underlying tight-binding band structure, i.e. the 



sensitivity to the choice of the Raman vertex in A\ g sym- 
metry, resolving the previously encountered problems in 
the symmetry analysis of the light scattering amplitude 
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On the basis of the observation of the collective SF 
mode in Y-123 and Bi-2212, we consider a bilayer model 
represented by a tight-binding band structure 

£k = — 2t(cos(fc;c) +cos(fc a )) +4t' cos(fc x ) cos(ky) — t±(k), 

with an inter-plane hopping given by M 



t±(k) = 2t±cos(k z )(cos(k x ) — cos(k y )Y 



(1) 



where k z is or tt for the bonding or anti-bonding bands 
of the bilayer, respectively. 

The spin susceptibility (xs) is modeled by extending 
the weak coupling form of a BCS superconductor in a 
d x 2_ y 2 pairing state to include antiferromagnetic spin 
fluctuations by a RPA form with an effective interaction 
U ; i.e. Xs(q, iw) = X°(q, (1 - [7x (q, iu)) where [§ 



X °(q,;a,) = iTr£ 



G(k, iu/)G(k + q, iJ 
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Tr denotes the trace and (3 = T 1 
Green's function in Nambu space 



G(k, iuj) is the BCS 
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with Ti (i=l,2,3) the Pauli matrices, fo the 2x2 unit ma- 
trix, £k = ek — [i and Ak = A [cos(/c 2; ) — cos(/c y )]/2. This 
form of the spin susceptibility contains a strong magnetic 
resonance peak at q = Q = (tt, tt, tt) which was proposed 
§ to explain the INS resonance at energies near 41 meV 
in Y-123 § and Bi-2212 @. Other forms for the spin 
susceptibility can be straightforwardly used within our 
model. However, the results are mainly determined by 
the collective mode at Q. Therefore we take the bilayer 
susceptibility for a representative calculation. 

The intensity of scattered light I(lu) in Raman exper- 
iments is proportional to the imaginary part of the re- 
sponse function for the effective density operator 
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in the long wavelength limit q — > 0. Specifically 
1(d) oc (1 + n(fi)) Imx(H + i0+) 

r-l/T 



,-tn 



MT T [p(r),p(0)]}, 
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with the Bose function n(u) and the time ordering oper- 
ator T T . 

The bare Raman vertex j(k) = T37(k) in different 
scattering geometries are classified according to the el- 
ements of the D 4h point group. For the limiting case of 
vanishingly small scattered (us) and incident (ui) pho- 
ton energies, it can be represented in the effective-mass 
approximation (EMA) 



(6) 



where e 1 and e s are the unit vectors for in-plane po- 
larizations (i.e. a,/3 G {x,y}) of the incoming and the 
scattered light, respectively. Using Eq. (||) and the bi- 
layer tight-binding dispersion Eq. [l] we obtain 



2i 7k (l + ^j- cos(^)[cos(fc x ) + cos(fc a )]^ (7) 



Big 

7 k = 

7k 19 = - 2t± cos(fc^)(cos(2fe x ) + cos(2fe s )) 

-4cos(fc 2; ) cos(fc„)[t' + 2tj. cos(fc z )] (8) 

where 7^' s = (cos(k x ) =p cos(fcj,))/2. However, the EMA 
has a questionable region of validity for all Raman mea- 
surements H on the cuprates since the incoming pho- 
tons have energy ~ 2eV, which is on the order of the 
bandwidth and the inter-band excitations according to 
local density calculations JTo| |. EMA based arguments 
in previous works about relative Raman intensities for 
different channels are therefore questionable ||. We 
hence consider other forms for the vertices as well which 
obey the proper symmetry transformations. For the 
A\ g geometry some symmetry compatible choices are 
7(k) = cos(fc x ) + cos(ky) and 7(k) = cos(fc x ) cos(k y ). 
These basis functions assign weight to different regions 
of the BZ and this is the reason why previous results 
for the Ai g response were particularly sensitive to the 
specific choice of the bare Raman coupling vertex. 

The spin fluctuations lead to an additional contribu- 
tion for the Raman response via a 2-magnon like process 
as shown diagrammatically in Fig. [I] pT| . Here, the SF 
propagator is incorporated in its RPA form for the bi- 
layer (as described above) by the ladder diagram series 
with an effective on-site Hubbard interaction U. 

We therefore write the Raman response function at 
finite temperature as the sum of a pair-breaking (PB) 
and a SF contribution 
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FIG. 1. a) Feynman diagrams for the Raman response 
function including pair-breaking quasiparticle excitations and 
"2-magnon processes". Dashed, wiggly, and solid lines repre- 
sent photons, spin fluctuations, and fermionic propagators. 7 
denotes the bare Raman vertex in a selected scattering geom- 
etry, b) Ladder series for the spin fluctuation propagator, c) 
Explicit diagram for the spin fluctuation contribution. 



with the Raman vertex specifying the scattering geome- 
try. In the limit q — > the diagram for the SF contribu- 
tion translates into 

= ~£ ^(q',^,^)x s (-q',-^) 

q' ,iuj 

x Xs (q', iu + in)V^(q', -iQ, -iu). (10) 

The vertex function V 7 (q', i£l, iu) includes the bare Ra- 
man vertex and is evaluated as 

7 7 (q',ffl,iw) = Trji ^ 7 (k)G(k,ia/ + ifi)f t7 

G(k + q, iu' + itt + iu)r UG(k, iu') 1 . (II) 



ifl, iu denote bosonic and iu' fermionic Matsubara fre- 
quencies. Similarly, x 7 ~f is evaluated as 



Tr 



(in) = — 7(k)G(k,iw')7(k)G(k,iw' + ifi). 

1 t ■ 1 
K.lttT 



(12) 

The total Raman response is calculated in the gauge 
invariant form which results from taking into account 
the long wavelength fluctuations of the order parameter 
to guarantee local charge conservation |§,|l2|. The total 
Raman susceptibility thus follows as 



X(ifi) = X77^) 



Xn(*^) 



(13) 



where X17 an d X11 are obtained by the replacement 
7(k) — > 1 in one or both bare Raman vertices in the 
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FIG. 2. "2-magnon process" contribution x'j-y i n the A± g 
geometry. The parameters used are defined in the text. 



vertex function Eq. ( |11| ) . The analytical continuation to 
the real axis is performed using Pade approximants jl3) . 

The band structure parameters are chosen for all 
the numerical calculations to be applicable to optimally 
doped systems: (n) = 0.85, t'/t = 0.45, t±/t = 0.1 §] 
while the gap has been chosen as Ao/t = 0.25. We 
have evaluated the Raman response at a temperature 
T/t = 0.08. 

Let's first consider xsf alone. In Fig. || we plot x'sf 
for the Ai g channel versus frequency, xsf is a convolu- 
tion of two spin susceptibilities (Eq. (|l0|) ) and is thus 
peaked at an energy near twice the magnetic resonance 
energy. An important point is that in the B\ g and £?2 9 
geometries the SF term introduces vanishingly small cor- 
rections to the total response, rendering the presence of 
the SF term important only in the K\ g geometry. This is 
due to the sharpness in momentum space of the resonance 
peak at Q in the SF propagator. In fact, if the transfer 
is taken only at Q, both the B\ g and Bi g contributions 
to xsf vanish identically, as can been seen from Eq. [II]. 
Therefore for these channels the response is given by the 
PB term alone. 

In Fig. H we illustrate the frequency dependence of 
%"(f2) in the Ai g geometry for three different values of 
the effective interaction with both the PB and SF con- 
tributions included. The shape of the Raman response is 
modified varying U/t and in particular the position of the 
resonance is shifted towards higher energies for increas- 
ing U. With the inclusion of the SF term we now obtain 
a peak slightly above Ao as observed experimentally in 
Y-123 and Bi-2212. 

A comment is in order on the relative magnitude of the 
SF and the PB term. Comparing Figs. || and [| it is clear 
that the SF term is much smaller than the PB term, but 
the effect of this new term is nevertheless visible since 
the backffow (the second term in Eq. (|l3|)) mixes in a 
non trivial way the two contributions. Since the SF term 
varies as U 4 in our model it starts to dominate for larger 
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FIG. 3. The total response in the A\ g channel. Solid, 
dashed and dotted lines correspond to U = 0, U/t — 1.3 
and U/t = 2.0 respectively. 



t7, leading to a shift in spectral weight out towards 2wr, 
which also changes with U. 

The experimental position of the peak is not sample 
dependent in the A\ g geometry as already mentioned, 
and is almost the same in different cuprates. On the 
other hand the theoretical description with the PB term 
alone is very sensitive to Raman vertex changes, which 
can produce variation of its position between Ao and 
2Ao not allowing for a comprehensive modeling of 
different cuprates. 

In Fig. ^ we address the problem of the sensitivity 
of the result to changes in the bare Raman vertex 7(k). 
To investigate the effect of changes of the vertex func- 
tion, we have calculated the final response using the three 
different forms for the vertex 7(k) = cos^) + cos(k y ), 
7(k) = cos(k x ) cos(k y ) and in the EMA which posses the 
correct transformation properties required by symmetry. 
In the first panel of Fig. 0, we have plotted the Raman 
response for [7 = 0, i.e. the PB contribution alone, and 
in the second panel for U/t = 1.3. All curves are renor- 
malized to their peak height to allow for an easier com- 
parison. Clearly the strong sensitivity to changes of the 
bare Raman vertex (first panel) is much reduced when 
the SF term is added (second panel). 

In Fig. ^| we compare the theoretical results with ex- 
perimental data on optimally doped Bi-2212 0]. Adding 
the SF contribution leads to a shift of the peak position 
from near ~ Ao for U = to higher frequencies, and 
thus to a better agreement with the experimental rela- 
tive peak positions in A\ g and B\ g geometries. For the 
fit we have adjusted t to achieve a good agreement with 
the B\ g channel, and then adjusted U to match the A\ g 
peak position. 

The value of t obtained from the fit is t — 130 meV. 
This value has to be compared with t « 105 meV, which 
results from the condition uir ks 40 meV. This slight 
discrepancy is most probably related to our simple mod- 
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FIG. 4. The total A\ g response for different vertices: 
7 obtained from effective mass approximation (solid line), 
cos(fc;c) + cos(fcy) (dashed line), cos(fc:r) cos(fe H ) (dotted line). 



eling of the propagators which neglects strong renormal- 
izations from interactions as well as impurities. 

From this work we conclude that including the SF con- 
tribution in the Raman response solves the previously un- 
explained sensitivity of the A\ g response to small changes 
in the Raman vertex. Also, within our model it is now 
possible to obtain the correct relative peak positions of 
the A\ g and the B\ g scattering geometry. Whereas the 
SF (two-magnon) contribution controls the A\ g peak, the 
Bi g and B2 g scattering geometries are essentially unaf- 
fected and determined by pair breaking processes alone. 
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